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Abstract 

The fluid-gravity correspondence provides us with explicit spacetime 
metrics that are holographically dual to (non-)relativistic nonlinear hydro- 
dynamics. The vacuum Einstein equations, in the presence of a Killing 
vector, possess solution-generating symmetries known as Ehler's transfor- 
mations. We apply these Ehler's transformations to vacuum metrics with 
hydrodynamic duals to develop a formalism for solution-generating trans- 
formations of Navier-Stokes fluids. Using this we provide examples of a 
linear energy scaling from RG flow under vanishing vorticity, and a set of 
Z2 symmetries for fixed viscosity. 
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1 Introduction 

In 1974, Damour [1], and later in 1986, Thorne et al. [2], considered an observer 
outside a black hole, interacting with (perturbing) the event horizon. Surprisingly, 
they found that the observer will experience the perturbations of the "stretched" 
horizon as modes of a viscous fluid possessing electric charge and conductivity. This 
inspired a host of works over the years but the connection between gravitational 
physics and fluids became sharper with the advent of the AdS / CFT correspondence 
when Policastro et. al [5] related the shear viscosity of M = 4 super Yang-Mills 
theory to the absorption of energy by a black brane. 

This was the start of using the holographic principle, that is the correspondence 
between gravitational theories on (d + l)-dimensional manifolds and <i-dimensional 
quantum field theories, as a tool for calculating hydrodynamic properties. 

More recently, there has been a set of works that directly relate solutions of 
Einstein's equations of a particular type to solutions to the Navier-Stokes equa- 
tions in one dimension less [6, 7, 8, 9, 10, 11, 12, 13, 14]. Later we will review 
the details of how this correspondence is derived but the essential flavour is as 
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follows. One writes down a very particular ansatz for the metric in d + 1 dimen- 
sions which has undetermined functions, v i (x, t), P(x, t) and parameter, v with the 
index i — 1, .., d — 1 i.e. over a (d— 1) subset of the d spacetime dimensions. Solv- 
ing the Einstein equations then constrains the functions v l (x, t), P(x, t) to give a 
set of second order nonlinear differential equations for v l (x, t), P(x, t). This set of 
equations are the Navier-Stokes equations describing a fluid in d dimensions with 
pressure, P(x,t), fluid velocity field v l (x,t) and viscosity v. Thus particular solu- 
tions to the Navier-Stokes equations provide particular solutions to the Einstein 
equations. 

Other recent work on the fluid/gravity correspondence may be found in [20, 
21, 22, 23, 24, 25, 26, 27, 28]. 

It has been known since Buchdahl [17] that for manifolds with isometries Ein- 
stein's equations have solution generating symmetries. That is there are "hidden" 
symmetries of the equations that map from one solution to the other. In fact there 
is a vast set of these as described by Ehler [18] and Geroch [19]. The question we 
wish to pose in this paper is whether the solution generating symmetries of Ein- 
stein's equations can lead to solution generating symmetries in the Navier-Stokes 
equations? 

The procedure to determine this will be as follows: 

• impose a Killing symmetry in a spacetime that admits the metric ansatz 
corresponding to the Navier-Stokes equations; 

• carry out Ehler's transformations that preserve the ansatz; 

• determine the induced transformation of the Navier-Stokes data i.e. v l , P, v. 

The procedure could immediately fail if there were no Ehler's transformations 
that preserve the metric ansatz required for the fluid/gravity correspondence. We 
will find that there are a finite set and we will be able to explore the transformations 
on the Navier-Stokes fields for different choices of Killing directions in spacetime. 

The paper will try to be as self contained as possible and so we begin with 
a review of the necessary ideas in fluids; the Navier-Stokes equation; the fluid 
gravity correspondence; and solution generating symmetries in general relativity. 
We will then carry out the procedure described above for spatial, timelike and 
null Killing vectors to see what solution generating symmetries they correspond to 
in the Navier-Stokes equation. We end with some comments and ideas for future 
work. A reader familiar with the formalism of hydrodynamics and the Navier- 
Stokes equation may wish to skip directly to section 3 where we carry out the 
duality transformation in the gravity dual to see the induced transformations. 
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2 Hydrodynamics 



The study of hydrodynamics is fundamental to vast areas of physics and engineer- 
ing, owing to its origin as the long-wavelength limit of any interacting field theory 
at finite temperature. Such a limit needs a consistent definition. Consider a quan- 
tum field theory where quanta interact with a characteristic length scale £ C orr ; the 
correlation length. The long-wavelength limit simply requires that fluctuations 
of the thermodynamic quantities of the system vary with a length scale L much 
greater than £ cmT , parameterized by the dimensionless Knudsen number 



For a physical description to be useful in non-equilibrium states, we naturally 
require that L remain small compared to the size of the system. This is usually 
satisfied trivially by considering systems of infinite size. 

The long-wavelength limit allows the definition of a particle as an element of 
the macroscopic fluid, infinitesimal with respect to the size of the system, yet con- 
taining a sufficiently large number of microscopic quanta. One mole contains an 
Avogadro's number of molecules, for example. Each particle defines a local patch 
of the fluid in thermal equilibrium, that is, thermodynamic quantities do not vary 
within the particle. Away from global equilibrium quantities vary between parti- 
cles as function of time r and spatial coordinates x, combined as x a = (r, x). The 
evolution of the quantities between particles is parameterized by a relativistic ve- 
locity u b (x a ), which refers to the velocity of the fluid at x a , rather than the velocity 
of an individual particle as it traverses a path through the fluid, parameterized 
by the coordinates. It is well known [29] that the thermodynamic quantities, such 
as the temperature T(x a ) and the density p{x a ), are determined by the value of 
any two of them, along with the equation of state. The evolution of the system 
is then specified by the equations of hydrodynamics in terms of a set of transport 
coefficients, whose values depend on the fluid in question. 

Fluid flow is in general relativistic in that the systems it describes are con- 
strained by local Lorentz invariance, and velocities may take any physical values 
below the speed of light. Applications at relativistic velocities are multitudinous: 
the dust clouds in galaxy and star formation; the flow of plasmas and gases in 
stars supporting fusion; the superfluid cores of neutron stars; the horizons of black 
holes are all described by hydrodynamics. Modelling black holes (and black branes 
in M/string theory) with hydrodynamics has now developed into a fundamen- 
tal correspondence of central importance to our present study, we elaborate in 
§1. Quark-gluon plasmas behave as nearly ideal fluids and are expected to have 
formed after the inflationary epoch of the big bang, are reproduced in collisions 
at the RHIC and LHC, and in the core of nuclear fusion reactors. Non-relativistic 
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fluids are equally ubiquitous, somewhat more familiar, and constitute an endless 
list of phenomena from the atmosphere to the oceans. 

2.1 Relativist ic fluids 

We begin with a discussion, adapted from [30], of the relativistic fluid described 
by the stress energy tensor T ab and a set of conserved currents Jf where I indexes 
the corresponding conserved charge. The dynamical equations of the <i-spacetime 
dimensional fluid are 

W a T ab = V a Jf = . (2) 

For an ideal fluid, with no dissipation, the energy-momentum tensor and currents 
may be expressed in a local rest frame in the form 

T ab = pu a u b + p{g ab + u a u b ) (3a) 
Jf = qi u a (3b) 

where p is the pressure and qi the conserved charges. The velocity is normalised 
to u a u a — — 1. The entropy current is given by (3b) with the charge q being given 
by the local entropy density. The conservation of the entropy current illustrates 
the non-dissipative nature intrinsic to zero entropy production. 

In a dissipative fluid, there are corrections to (3). First, we must first take 
into account the interrelation between mass and energy to define the velocity field 
more rigorously. This is achieved by using the Landau gauge, which requires that 
the velocity be an eigenvector of the stress-energy tensor with eigenvalue the local 
energy density of the fluid (this is satisfied by the velocity normalisation for the 
ideal fluid). If the stress energy tensor gains a dissipative term IP 6 , and the current 
a term T a , this reads 

W b u a = T a u a = 0. 

Dissipative corrections to the stress tensor are constructed in a derivative expansion 
of the velocity field and thermodynamic variables, where derivatives implicitly 
scale with the infinitesimal Knudsen number (1). Recalling that the equations of 
motion for the ideal fluid are composed of relations between these gradients, we 
may express IT afc purely in terms of the derivative of the velocity. This can be 
iterated to all orders in the expansion. Now, the derivative of the velocity may be 
decomposed using the acceleration A a , divergence 9, a symmetric traceless shear 
a ab , and the antisymmetric vorticity u ab into the form 

VV = —A a u b + a ab + u ab + — l —9P ab , 

a — 1 
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and P ab = g ab +u a u b is a projection operator onto spatial directions. In the Landau 
frame, only the divergence and shear can be built purely from velocity gradients, 
and so define the leading order contribution to the stress-energy tensor. A similar 
analysis for the charge current retains the acceleration, and if one includes the 
parity-violating pseudo-vector contribution 

t = e bcd a u b V c u\ 

the leading order dissipative equations of motion for a relativistic fluid are (2) with 



T ab = pu a u b + pP ab - 2r]a ab - (6P ab (4a) 
Jf = qiU a - XuP ab V b qj - Qj£ a - nP ab V b T, (4b) 

where rj and ( are the shear 1 and bulk viscosities respectively, (u is the matrix 
of charge diffusion coefficients, 7/ indicates the contribution of the energy density, 
and 0/ the pseudo-vector transport coefficients. The transport coefficients have 
been calculated in the weakly coupled QFT in perturbation theory, whereas in the 
strongly coupled theory, a dual holographic description may be employed, see e.g. 
[3]. 

Returning to the discussion of the entropy current, we now expect that the 
entropy will increase in fluid evolution to equilibrium, consistent with the second 
law of thermodynamics. This translates to the requirement that 

W a J a s = ^a ab a ab > 0, 

which is satisfied iff the shear viscosity is positive, a condition satisfied for a large 
class of QFTs [4]. 

1 We use the traditional notation of [4, 12, 29] rather than that of [8, 9]. We define the 
kinematical viscosity v in §2.2. 
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2.2 Non-relativistic Navier- Stokes fluids 

We now take the non-relativistic limit, denned by long distances, long times and 
low velocity and pressure amplitudes, of the relativistic fluid on a spacetime lightly 
perturbed in the following manner [12]. Consider an arbitrary perturbation 

Gab = (jab + H ao 

of the background metric 

g a b = -dr 2 + g ij dx l dx J 
(this form is purely a choice of coordinate system), where 

H TT = e 2 h TT (e 2 r, ex 1 ) 

H Ti = eai(e 2 T,ex l ) (5) 
Hij = e 2 hij(e 2 T,ex l ). 

Treating the fluid on G a b as an effective forced flow of that on g ao , whilst preserving 
the normalisation u a u a = — 1, the velocity 

u a = -. 1 . =(1,V% 

density and pressure may be expanded about an equilibrium configuration (p , V$, po) 
of a stationary, uniform fluid at rest, as 

V i = e^(e 2 r, ex*) 

(6) 



Po +Po 



e p{e r, ex 1 



with an infinitesimal scaling parameter e. The scaling of the density is fixed by the 
equation of state. The scaling of time r and space x % are implicit. In this limit, 
the temporal component of the relativistic Navier-Stokes equations becomes 

V a T aT = e 2 (p + Po)V^ + 0(e A ) (7) 

where Vj is the covariant derivative with respect to the spatial metric Sepa- 
rating the gauge field di into its pure curl and divergence parts 

hi = a,i + ViX, 

such that Via 1 = 0, and defining an effective pressure 

P = v --h +^ 

6 



the spatial components yield 

W a T ai =e 3 (p + Po) [V,P + d T Vt + v j V jVi - u(V j V jVi + R ijV j ) 
+ d T a t + vJf Jt ] + (D(e 5 ) 

where 

v = V/ (Po + Po) 

is the kinematical viscosity, fy = didj — djdi is the field strength of aj (equivalently 
Sj) and Rij is the Ricci tensor on If a fluid carries conserved charges, the 
equations V a Jf = will also yield incompressibility (7). Taking the hydrodynamic 
limit e — > 0, these give the incompressible non-relativistic Navier-Stokes equations 



with a forcing function due to an external electromagnetic field, 

d T v i — v[d 2 Vi + RijV j ) + diP + v j djVi = —d T cii — v j fji (9a) 

dtf = 0. (9b) 

Ideal fluids are described by Euler's equations, obtained by setting v = 0. We 



will mostly be concerned with fluid flow on flat space (Rij = 0) in the absence of 
external forces, where H a b, thus a iy are zero. 

2.3 The Navier-Stokes fluid on a Rindler boundary 

A metric dual to the non-relativistic incompressible Navier-Stokes equations was 
first developed in [8] on the Rindler wedge. This included terms up to third order 
in the e expansion as in (6). An algorithm for generalising this metric to all orders 
was subsequently developed in [9], though terms calculated beyond third order are 
not universal. They receive corrections from quadratic curvature in Gauss-Bonnet 
gravity [11]. We summarise the construction in [9] here. 
Consider the surface S c with induced metric 

-f ab dx a dx b = -r c dr 2 + dxidx* (10) 

where the parameter y/r c is an arbitrary constant. One metric embedding this 
surface is 

ds 2 = -rdr 2 + 2drdr + dx { dx\ (11) 

which describes flat space (fig. 1) in ingoing Rindler coordinates x M = (r, x l ,r), 
defined in terms of the Cartesian chart (t, x\ z) by 

z 2 -t 2 = 4r z + t = e T/2 . (12) 

The hypersurface S c is defined by r = r c where r is the coordinate into the 
bulk. Allowing for a family of equilibrium configurations, consider diffeomorphisms 
satisfying the three conditions 
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Figure 1: The past H~ and 
future 1-L + horizons define 
the boundary of the Rindler 
wedge. Grey lines demon- 
strate lines of constant r 
(curved) and r (straight). 
Long-wavelength perturba- 
tions of the hypersurface S c 
are described by the equations 
of hydrodynamics. 



i) The induced metric on the hypersurface S c takes the form (10) 

ii) The stress tensor on S c describes a perfect fluid 

iii) Diffeomorphisms return metrics stationary and homogeneous in (r, x l ). 

The allowed set is reduced to the following boost, shift and rescaling of x M . First, 
a constant boost 



V^r -»- j(^F c r - fax*), x l x l - + (7 - l)^r^, (13) 

where 7 = (1 — /3 2 ) -1 / 2 and = r c ~ 1 / 2 v j. Second, a shift in r and a rescaling of r, 

r^r-r h , r ^ (1 - r h /r c y 1/2 r. (14) 
These yield the flat space metric in rather complicated coordinates, 

ds 2 = — ; — I v 2 1—) = =drdr , 2 4 =dxMr 



l-v 2 /r c V l-r h /rj ^1 - r A /r c r c - V r c 

+ a» / r-M di>dr + / _ ^ / r-*J) \ ^ 

1 - v 2 /r c \r c - r h J \ r 2 c {l - v 2 /r c ) \ 1 - r h /r c J J 

(15) 

The Brown- York stress tensor on S c (in units where 167rG = 1) is given by 

T ab = 2{K lab -K ab ), (16) 
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where 

K ab = \{C nl ) ab , K = K\, 

are the extrinsic curvature and its mean, and is the spacelike unit normal to 
the hypersurface. 

By imposing that the Brown- York stress tensor on S c gives that of the stress- 
energy tensor of a fluid we can identify the parameters of the metric (15) with the 
density, p, pressure, P and four-velocity u a of a fluid, as follows: 

P = 0, p= * , u a = -=L=(l,Vi). (17) 

V r c - r h \Jr c - v 2 



The Hamiltonian constraint 



R^rfn" = 



on E c yields a constraint on the Brown- York stress tensor 

dT ab T ab = (T a a ) 2 . 

When this constraint is applied to the equilibrium configurations described above, 
one finds the equation of state is p = (as above), or p = —2d(d — l)p which 
occurs for a fluid on the Taub geometry [14]. 

Promoting V{ and p to slowly varying functions of the the coordinates x a , and 
regarding Vi(r,Xj) and p = Tc 1 ^ 2 + Tc 3 ^ 2 P(t, Xj) as small perturbations about 
equilibrium yields the metric 

ds 2 = — rdr 2 + 2drdr + dx l dxi 

2(1--] UfdarMr - 2— darMr 



r 

+ U-- 

r r 



( v 2 + 2P)dT 2 + V -^dx i dx> ' l ' 2 + 2P 



• ( - I " ) drdr + 0(e 3 ) 

(18) 



which satisfies the Einstein's equations to 0(e 3 ) if u» satisfies incompressibility, 
d{V % = 0(e 3 ). Corrections appear in powers of e 2 , so this is the complete metric to 
second order. 

The metric may now be built up order by order in the hydrodynamic scaling. 
Assume one has the metric at order e n , where the first non- vanishing component 
R$ of the Ricci tensor appears at order n. By adding a correction term g$ to 
the metric at order n, resulting in a shift in the Ricci tensor 8R$ , and requiring 

R$+6R$ = 0, (19) 
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the vanishing of the Ricci tensor is guaranteed to order n. Recalling that, in the 
hydrodynamic scaling, derivatives scale thus, 

<9 r ~e°, di~e\ <9 r ~e 2 , (20) 

one sees that corrections 5R$ at order n will appear only as r derivatives of gffi. 
It is shown in [9] that, using the Bianchi identity and the Gauss-Codacci relations, 
integrability of the set of differential equations (19) defining SRpJ in terms of g$ 
is given by imposing the momentum constraint, equivalent to the conservation of 
the stress tensor on S c , 

RaX = V a T ab |g = 0, (21) 

which is precisely the fluid equations of motion, to order n. 

The perturbation scheme contains several degrees of freedom. The gauge free- 
dom of the infinitesimal perturbations 

9$ 9$ + d&P + duVt* 

for some arbitrary vector ip^ n '(r,x,r) at order e n , which may be fixed by de- 
manding that g rfl is that of the seed metric to all orders in e. The a: a -dependent 
functions of integration from equation (19) may be fixed by imposing the bound- 
ary form (10) of the metric on E c , and also requiring regularity of the metric at 
r = 0, which in this construction translates to the absence of logarithmic terms in 
r. Corrections to the bulk metric under these conditions then become 

g£) =(1 - r/r c )F^{r,x) + H dr' f° dr"(R% ] - rR%> - 

J r J r f 

g %> =(1 - r/r c )F} n) (r, x) - 2 H dr' H dr" R^ (22) 

Jr Jr' 

pr c i pr' 
An) _ o A J AJI E>(«) 



C=-2jf dr^jf dr»R%>, 

where the Fa ( T , %) comprise of the remaining integration functions, and the final 
degree of freedom; field redefinitions of Svj- and 5P^ at order e n . F^(t,x) is 
related to redefinitions of the fluid velocity and is fixed by the isotropic gauge 
condition Pj^Tb c u c = 0. Fr (r,x) is related to redefinitions of the pressure and is 
fixed by defining the isotropic part of to be 

lT rOPlC = (4= + ( 23 ) 
\ V c r c 

to all orders. 
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Applying the perturbation scheme to the seed metric yields to third order, 



ds 2 



rdr 

1 - 

fr 2 



2<\rdr + dx'dxi - 2 | I 
v 2 + 2P)dr 2 -} 



r 
r. 



2 ' "^diW 



WjdxMr — 2— dxMr 
+ | ) drdr 



r,. 



(24) 



1 - 



v 2 + 2P 



r c V W V r c 

which satisfies the vacuum Einstein equations if 

r r 3/2 V a T a J Er = d T Vi - r c d 2 Vi + diP + vWm 



which are the Navier-Stokes equations with kinematical viscosity 



v 



(25) 



(26) 



The corresponding corrections to the Navier-Stokes equations follow from con- 
servation of the stress tensor on E c . Vector and scalar quantities are odd and 
even orders respectively in the scaling e. Accordingly, corrections to the scalar 
incompressibility equation appear at even orders, and to the vector Navier-Stokes 
equations at odd orders. 

It is shown in [8] that the hydrodynamic limit is equivalently described on 
the Rindler wedge by taking the boundary to infinity r c — > oo, which effectively 
renders deviations from the horizon infinitesimal. 



2.4 The Navier-Stokes fluid on a Schwarzschild sphere 

In [10], Bredberg et al. construct the (d+ l)-dimensional metric dual to the non- 
relativistic fluid on a <i-dimensional hypersphere (S denotes Schwarzschild) a 
distance r s > from the horizon of a Schwarzschild black hole. The Dirichlet 
boundary conditions used on the Rindler space are not consistent here. On Ej;, 
one instead demands conformal flatness where the conformal factor is allowed to 
fluctuate, and constant mean curvature K defined with the unit tangent orthogonal 
to Eg. The perturbative expansion is parameterized by the inverse square of the 
curvature, 

-2 



r. 



1/AK 2 



in the limit K — > oo. The metric is perturbed from the Schwarzschild metric 



ds 2 

uo Schwarzschild 



2m \ 
r J 



(d-2) - 



dt 2 + 



2m 
r 



(d-2)' 



df 2 + r 2 d^_!. 



(27) 
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where dfl 2 l _ 1 is the metric on the unit (d — l)-sphere. In ingoing null coordinates 
(r, r) defined by 

(d-2)a( f ? dr' \ 2m ^ n . 

t = K - ' \t+ 7~, — ) j r = — — — (f-2m), 28 

2m y J i_(2rny d -^J a(d-2) y h y ' 

the leading terms in the perturbed metric obeying the boundary conditions are 

j2 r l2 (d — l)(d — 2)r(5 + r) 2 nJ , A 2 ir ^ 2 

ds 2 = dr 2 + ^dr 2 + 2drdr + 4m 2 dfiL 1 

r s 8m 2 

+ (1 - r)[v 2 dr 2 - 2v i drdx i ] - 2rPdr 2 + r s (2(d - 2)r + 8m 2 P)dfi^ 1 

+ (1 - r)v i v j dx i dx j - (r 2 - 1)(V 2 ^ + (d - 2)R ji v j )drdx i - 2^drdx i 

+ (v 2 + 2P)drdr + 2(1 - r^drd^] + r 2 ^XidxMr + A^dxMx^ + 0(rf ) 

(29) 

for functions ((j>i(x a ),Xi(r,x a ), \ij(r,x a )). The Ricci tensor on the (d — l)-sphere 
takes the form i?^- = r~ 2 (d — 2)gij. The perturbed geometry satisfies Einstein's 
equations to O{r s ) if (v iy P) satisfy the incompressible Navier-Stokes equations on 
the sphere with viscosity v — 1. In addition, the divergence of <pi(x a ) is fixed in 
terms of (v i: P), and (2V l x« — <9 P A-) is determined in terms of some function regular 
in r. 

As Bredberg et al. point out, this spacetime points to a connection between 
famous conjectures on either side of the correspondence. The cosmic censorship 
conjecture states that singularities do not form on or outside event horizons in a 
spacetime arising from smooth initial data, while global existence for the Navier- 
Stokes similarly suggests that singularities do not form from smooth initial velocity 
fields. One can see that if a singularity were to form from smooth initial data, 
in the fluid velocity on S c outside the horizon, this would simultaneously violate 
both conjectures. 

3 Duality in the context of holography 

The defining equations in general relativity are the Einstein field equations, and in 
the non-relativistic limit of hydrodynamics, the Navier-Stokes equations (9). Each 
is a set of non-linear partial differential equations describing the evolution of two 
fundamentally inextricable systems, in themselves exhibiting fantastically varied 
phenomenology. When approaching any complex physical system with a view to 
finding solutions, it is often advantageous to consider the symmetries, intensively 
studied in both of these systems since their conceptions. Beyond diffeomorphisms, 
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the search in gravity has in general been somewhat limited [15, 36], however in the 
presence of a spacetime isometry, the symmetry group becomes remarkably large 
[35], particularly for vacuum spacetimes. For symmetries of the Navier-Stokes 
equations see [31], and with regards to the conformal group [12, 32]. In the light 
of the fluid/gravity correspondence, one may ask whether the symmetries of these 
systems are linked. In [37, 38, 39], they apply known symmetries of the Einstein 
equations to spacetimes with perfect fluid sources, constructing new spacetimes 
with the same equation of state, though not within a holographic framework. By 
drawing on the tools provided by these works and those in holography, we hope 
to develop a more general approach to the problem. 

The bulk provides an additional valuable degree of freedom, where the bound- 
ary sets the scene for the fluid evolution on the induced geometry. Moreover, we 
are now free to exploit the symmetries of the more extensive yet simpler vacuum 
geometries. It is these symmetries which we intend to holographically project to 
the fluid. In particular, we are interested in dualities between solutions to the 
Navier-Stokes equations arising from dualities between solutions to the vacuum 
Einstein equations: dual metrics yield dual fluid configurations. 

In this section, we discuss the work leading up to the spacetime Ehler's sym- 
metry group of the vacuum field equations with zero cosmo logical constant. We 
continue in §4 to apply this to fluids on the boundary of the Rindler space, thus 
deriving solution generating transformations of the fluid velocity, pressure and vis- 
cosity (the latter defining the RG flow). We offer in §4.1.1 a selection of example 
transformations including RG flow for zero vorticity fluids (where one may relax 
this constraint), and Z2 transformations for fixed viscosity which we show in §4.2.2 
in fact lie outside the spacetime Ehler's group. 

3.1 Symmetry groups of the Einstein equations 

Understanding the properties of the Einstein field equations has long been a subject 
of great theoretical interest, a sensible starting point being the inherent symme- 
tries involved. To this end, Buchdahl [17] derived a form of duality in vacuum 
spacetime metrics, where an n-dimensional vacuum metric stationary with respect 
to a coordinate x s : 

9ixu, s = 0, g^ s = /I, v G {0, . . . , n} 
generates a dual vacuum metric 

V = ((9ss) 2/{n ~ 3) g a b, (^r 1 ), a, b e {fi ^ s}. 

It is this solution generating property of spacetime isometries we wish to apply 
to solutions of the Einstein equations and holographically map to hydrodynamics. 
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We have, however, a considerably larger symmetry group at our disposal. The 
works of [16, 18, 19] develop the concept culminating in the generalized Ehlers 
symmetry group of the Einstein equations also for non-vacuum spacetimes. An 
extension exists [33, 34] to dualities between vacuum spacetimes and those with 
electromagnetic backgrounds described by the Einstein-Maxwell equations, of rel- 
evance for magnetohydrodynamics. We are concerned with the vacuum subgroup 
of the generalized Ehlers group, which includes the spacetime Ehlers group. 

3.2 The spacetime Ehlers group 

Let us begin by defining the Killing vector. Consider a manifold A4, endowed 
with a metric g. We define the set of Killing vectors £ G TM. as the set of 
diffeomorphisms preserving the metric. This is equivalent to the statement that 
that Lie derivative of the metric along £ = £ M <9 M vanishes: 

OCtfW = i p 9^, P + g„p + e,u g w = o. (30) 

Consider a four-dimensional Lorentzian spacetime manifold with metric g = 
(g^u) satisfying the vacuum Einstein equations, and admitting some Killing vector 
£. The twist potential 

= V-det(p) e^'VT (31) 

and Killing vector norm 

A = -e^ (32) 

give the Ernst one-form 

o"^ = = V M A - iujf, (33) 

for some scalar q (exactness is guaranteed by vanishing Ricci tensor, see [40] p. 164). 
Define a self-dual two form 

= (1 + 1*) V M £, , (34) 

where * is the Hodge dual operator. Define a one form W satisfying 

V {I1 W U] = -2 7 »[(7 ? + %S)^\ (35a) 
CW a + 1 = (i^ + 5)(-i^+5) = ( 7 A) 2 + (ju + 5f (35b) 

where a bar denotes complex conjugation, and 7 and 5 are non-simultaneously 
vanishing real constants, which as a pair fix the gauge of W. There exists a family 
of metrics, also with Killing vector £, which satisfy the vacuum Einstein equations, 
given by the map 

9r -> M£> W> 9) = ^ V - 2^W U) - ^ 2 W,W V (36) 
where Q 2 = £ a W a + 1 > 1 and the inequality holds over the whole geometry. 
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4 Solution-generating the Navier-Stokes fluid 



It is clear that the dual metric h(£, W, g) will in general take a different functional 
form from g. Consider the subset of transformations (36) preserving the functional 
form of g((p) up to redefinitions <f) — > <fi of some set of parameters. Vanishing of 
the Ricci tensor will yield a relation common to both metrics, in the parameters 
(f> and their redefinitions. 

In the case of the Rindler metric dual to the incompressible Navier-Stokes fluid, 
define these parameters by the fluid velocity Vi, pressure P, and boundary position 
r c within the bulk. The dual metric is parametrised by a set of dual parameters 
Vi, P and f c , denoted by On satisfying the vacuum Einstein equations on E c , 
now at r = f c in the dual geometry, the dual metric will yield the incompressible 
Navier-Stokes equations in the dual parameters 

diVi = (37a) 
d T Vi + diP + v k d k Vi - f c d 2 Vi = 0. (37b) 

Vitally, if (v^ P) satisfy the Navier-Stokes equations with viscosity v = r c , then 
the dual velocity and pressure (di, P) represent a dual set of solutions for viscosity 
v = r c . That is, a given subset of the spacetime Ehlers transformation acting 
on the fluid metric (24) obeying some Killing isometry corresponds to solution- 
generating transformations of the velocity and pressure, and RG flow parametrised 
by r c , of an incompressible Navier-Stokes fluid. 

At this point, we may generalise somewhat. The Rindler metric is just one 
fluid metric supporting flat background geometries on the boundary. We therefore 
wish instead to retain the common features of such metrics; the metric gauge 
g^ r , and the flat boundary metric which, by transforming to a suitable coordinate 
system, may always be expressed in the form (10). There is of course the additional 
freedom with vacuum metrics to introduce a constant conformal metric rescaling, 
varying under the duality, which we choose to neglect in our calculations. The 
equation we wish to solve is thus 



g^u(r c , Vi, P) ->■ h^it, W, g) = g^ = g^(r c , v h P). (38) 



where 



v(x a ) 2 + 2P(x a ) „ Vi(x a ) . . , 

9rr = 1 + V \, — , # r = — g rr = 39a 

2r c r c 

9ab\f c = lab, where 7 TT = -f c , % b = >y ab a ^ t ^ b. (39b) 
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4.1 Transforming the fluid 

We are provided in (39) with sufficient information to derive the fluid transforma- 
tion via the form of the one-form W^. Preserving the vanishing of the g rr = g rr = 
component of the metric we find, directly from (36), the two possibilities 

W r = -2a£ r tt 2 /X where a = 0,1. (40) 

One may obtain an expression for W a by contraction of (36) with the boundary 
indices (a, b, . . .) of the Killing vector 

X/Q 2 + (1 - 2a)£% A • 1 } 

This expression is uniquely defined only at the dual boundary E c , where we have 
defined the form of g a b and W a becomes independent of the dual fluid velocity 
and pressure. These expressions diverge for null Killing vectors, where A = 0. We 
cover this case shortly. 

One can see that the definition of the fluid is given, to all orders in e, by g ar = 
9ar\r c - Consequently, the fluid duality will be given by transformation of these 
components. Evaluation at f c is necessary in order to circumvent the ambiguity 
in the dual metric, and also provides explicit fluid transformations. We begin 
with Killing vectors null at the dual boundary, A|r c = 0, where one finds from 
contraction of (36) with the Killing vector, which yields 

^W, = C(g,u-g,u), (42) 



the following transformation, 

'€ b (.9ab - lab) 



9ar 9ar ~^ 



(43) 



accompanied by the preservation of a null Killing vector, ^^ u g^ v 
For non-null Killing vectors we employ the relations 



C{9ar - (1 - 2a)g ar ) = 0, (44) 
derived by comparing (40) and (42), and 

x/n 2 = -eCg»v, (45) 

found from contraction of (36) twice with the Killing vector. Inserting W r (40) 
and W a (41) into the Ehlers transformation (36) and employing (44) and (45), one 
finds 

-Xg ar + U^-2a)e% b -ia) 



9a 



ee rf 7cd + (l-2a)e^ 



(46) 

fc 
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4.1.1 Energy scaling invariance from a bulk isometry 

We begin with an example a (null) Killing vector into the bulk, £ = ^ r {x^)d r . 
Integrability of the Killing equations requires 

2vjd[iVj] = -r c d 2 Vi (47) 

along with vanishing vorticity, which with incompressibility implies (47). Trans- 
formation (43) gives 

* = P = f -l P+ t,( 1 . f A^ (48) 
r c r c r c \ r c J 2 

It is trivial to show that the pair (di, P) satisfy the incompressible Navier-Stokes 
equations if (vi, P) do so for velocities satisfying (47) alone - vanishing vorticity 
imposes unnecessary constraint and removes the dissipative term from the fluid 
equations. 

It is interesting to consider the problems of existence, uniqueness and regularity 
of the Navier-Stokes in this case. The divergence of (47) yields a vanishing mean 
square vorticity which ensures the class of solutions (v^ P) generated by (48) are 
regular. With respect to existence, the kinetic energy scales by a factor r c /r c and 
thus is bounded if there exists any solution satisfying (47) where the energy is 
finite. 

We note that the twist vanishes, and there exist values of 7/5 for which Q 2 = 
and the dual metric exhibits curvature singularities. We do not know how, if at 
all, this translates to the dual fluid. 

4.1.2 The timelike Killing vector 

One might expect from the timelike Killing vector £ = d T a duality of the form 

Vi = -Vi, P = P, f e = -r c (49) 

enacting time-reversal of the fluid but this is not the case. This is explained by 
noting that time-reversal is enacted by redefining the viscosity by v = ±r c [8] 
rather than by changing r c itself. This is because sending r c — > —r c brings the 
fluid outside the causal region of the spacetime. 

4.2 Fixed viscosity dualities 

We turn to fixed boundary (viscosity) transformations, where r c = r c . For Killing 
vectors null at the dual boundary, or if a = 0, one recovers the identity. For 
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non-null Killing vectors with a — 1, one finds 

9ar ~ 2U b lab/^ f lef 



9ar 

which constitutes a Z 2 group. 



1 - ezr/^ d ic d 



i 



(50) 



4.2.1 Spacelike Killing vectors 

Consider a generic space-like Killing vector £ = i k d k . Under (50), the pressure is 
preserved, while the velocity transforms as 

v i = Vi -2t£^, (51) 

which is a reflection in the hyperplane normal to the Killing vector and containing 
the point at which the velocity is defined. 

Translational isometry 

Consider £ = c k d k where c k are normalised to ^ k c 2 k = 1, and the corresponding 
isometries are c k d k Vi = c^d^P = 0. The dual fields are 

Vi = Vi- 2ciC k v k P = P. (52) 
The incompressibility condition 

diVi = diVi - 2cidiC k v k = 0, (53) 

and Navier-Stokes equations 

d T Vi + diP + v k d k vi - r c d 2 Vi =(5 ik - 2ciC k )(d T v k + d k P + v J djV k - r c d 2 v k ) 

+ 2ciC k d k P - 2c j v j c k d k (v i - 2ciCiv t ) = 0, 

are satisfied by the incompressible Navier-Stokes equations in the original fluid pa- 
rameters along with the isometries. This is valid for fluids of arbitrary dimension. 

Rotational isometry 

Consider a Killing vector £ = — x 2 d\ +Xi<9 2 corresponding to a rotational isometry 
in a rf-dimensional fluid. In polar coordinates 

xi = pcos9, X2 = psin9 x k > = x k >, (55) 
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where primed (') indices run from 3 to (d — 1), and the Killing vector becomes 
£ = dg, the isometries are 

devi = -v 2 d e v 2 = vi d e v k , = d d P = 0. (56) 

Solutions satisfy 

Vi = xifj,(T,x k >,p) - x 2 rj(T,x k ,,p) v 2 =x 2 fx(r,x k/ ,p)+x 1 r](T,x k/ ,p) (57) 
with dual velocities 

vi = x 1 p(r,x k >,p) + x 2 r)(r,x k/ ,p) v 2 = x 2 p(r, x k >, p) - x^t, x k >, p), (58) 
that is, the transformation sends rj — >■ —rj (equivalently 9 — > —9). The incompress- 



ible Navier-Stokes equations for the original fluid may be expressed as 

= (2 + pd p )p + d k ,v k , (59a) 

= d T p + p~ l d p P + /i 2 - i] 2 + ppdpp + v k >d k rp (59b) 

-r^p^dpp + d^ + d^d^p) 

= d T r] + 2pr] + ppd p r] + v k >d k <r] - r c (3p _1 9 p ?7 + djr? + d k ' dyr]) (59c) 

= <9 T v fe / + d k >P + ppd p v k , + VfdfV k > - r c p~ l dp(pd p v k >) (59d) 



It is clear from the parity of these equations in r] that if there exists a fluid solution 
defined in terms of a pair (p, rj) by (57), then there also exists a solution param- 
eterized by the pair (p,—rj). That is, the dual fluid satisfies the incompressible 
Navier-Stokes equations. Again, this is valid for fluids of arbitrary dimension. 



We provide an example with the three-dimensional fluid solution 

vi — A ( Xl - x 2 e~ 2A{T ~ TQ) ) v 2 = A (x 2 + x ie - 2A(T - To) ) (60a) 

v 3 = B exp (aA{t - r ) + ^pj - e 2Ar J dr' q{r')e- 2M ' - 2Ax 3 (60b) 

P = ]-A 2 p 2 (e~ AA{T - To) - 1) + q(r)x 3 - 2A 2 xj (60c) 



which satisfies the isometries (56). Here, A, B and r are arbitrary non- vanishing 
constants and q{r) is an arbitrary function of time. The duality is equivalent to 
sending r — > r + in/2A. 

4.2.2 Group structure for fixed boundary 

Consider the fixed boundary transformations in the case £ r = 0. The constraints 
(35a) on d[ r W a ] are dependent on radial derivatives of the dual metric, which we 
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have not denned, and so remain arbitrary. The RHS of (35a) is closed [16], so 
it remains to constrain 0r o W&i| rc , which vanishes. Now, we may write the Killing 
vector norm in the form 

- A = 7 a« b + (r c - r){f (flW + 0(e 4 )) + ^iffkr + 0(e 3 ))} 2 , (61) 

which by (44) is preserved, up to fourth order terms, under the duality. Equation 
(45) then yields Q 2 = l + (9(e 4 ), and the constraints are satisfied for 7 = 0, 5 = ±1. 

One begins in four dimensions [16] with the spacetime Ehlers group isomor- 
phic to SX(2,R), parametrised in terms of the Mobius map of the Ernst scalar 
q —7- + z7')/(i7<? + 5) by the real constants (7, 5, 7', 5') satisfying Y7 + 5' 5 = 1. 
For fixed boundary, 7 = and 5 = ±1 in the Mobius map result in a constant 
imaginary shift in the Ernst scalar which can be absorbed as an arbitrary inte- 
gration constant of the twist. In this way, the metric transformation is decoupled 
from the SL(2,M.) action. 

Transformation (36) forms a group acting on metrics in arbitrary spacetime 
dimension [16], which is not in general vacuum-to- vacuum, but contains the four 
dimensional spacetime Ehlers group as a particular case. The fixed boundary fluid 
transformations arise from vacuum-to- vacuum maps within the generalized Ehler's 
group, but outside the spacetime Ehler's group. For those examples in §4.2.1, ex- 
plicit calculation of the action of (51) on the Ernst scalar in four spacetime dimen- 
sions yields q — y q(l + 0(e A )) (provided the arbitrary integration constant of the 
twist vanishes), which provides some hint as to the form of this spacetime trans- 
formation. Note, in the presence of non-vanishing twist, this complex conjugation 
is not found in the Mobius map for any values of the group parameters. 

5 Discussion 

We have demonstrated how solution-generating transformations of the vacuum 
Einstein equations in the presence of a Killing vector may be applied to spacetimes 
holographically dual to hydrodynamics. Our focus has been on the incompressible 
Navier-Stokes fluid dual to vacuum Rindler spacetime, where we have uncovered 
a selection of dualities: a linear energy scaling symmetry for solutions with van- 
ishing vorticity (this constraint may be relaxed to (48)), deriving from RG flow 
of the fluid hypersurface through the bulk, and a Z 2 group of transformations for 
fixed viscosity (boundary) with explicit examples of reflection-like symmetry in 
translational and rotational fluid isometries. These transformations may not be 
remarkable from the perspective of hydrodynamics but it shows how part of the 
Ehler's transformations can survive holography and give rise to transformations 
in the fluid dual. One can then be inspired to try the Harrison transformation 
which is a solution generating transformation in Einstein-Maxwell theory to give 
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new transformations in magnetohydrodynamics. A holographic relation of the sort 
described here has been constructed for magnetohydrodynamics in [41, 42]. The 
Harrison transformation in the bulk may then lead to nontrivial transformations 
in magnetohydrodynamics transforming between fluid velocity and magnetic po- 
tentials; this is the subject of current work. One can also study the dimensional 
dependence of these solution generating transformations. In recent work [43], the 
difference in scaling for turbulence between three and four dimensions was studied 
holographically with a large difference in qualitative behaviour. One could also 
examine backgrounds associated with nontrivial chemical potential in the dual; for 
example rotating or charged black hole backgrounds. 

Essentially, in this paper we wish to open up the use of gravitational solu- 
tion generating symmetries in holography. It is encouraging that this did not 
immediately fail and one could preserve the fluid metric ansatz with some residual 
transformations surviving, yet it is intriguing that these transformations did not 
give anything particularly new. The results for magnetohydrodynamics may prove 
more significant. 
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